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ABSTRACT 

The radii of giant planets, as measured from transit observations, may vary with wavelength 
due to Rayleigh scattering or variations in opacity. Such an effect is predicted to be large 
enough to detect using ground-based observations at multiple wavelengths. We present de- 
focussed photometry of a transit in the HAT-P-5 system, obtained simultaneously through 
Stromgren u, Gunn g and r, and Johnson / filters. Two more transit events were observed 
through a Gunn r filter. We detect a substantially larger planetary radius in u, but the effect 
is greater than predicted using theoretical model atmospheres of gaseous planets. This phe- 
nomenon is most likely to be due to systematic errors present in the u-band photometry, stem- 
ming from variations in the transparency of Earth's atmosphere at these short wavelengths. We 
use our data to calculate an improved orbital ephemeris and to refine the measured physical 
properties of the system. The planet HAT-P-5 b has a mass of 1.06 ± 0.11 ± 0.01 Mjup and a 
radius of 1.252 ± 0.042 ± 0.008 Rjup (statistical and systematic errors respectively), making 
it slightly larger than expected according to standard models of coreless gas-giant planets. 
Its equilibrium temperature of 1517 ± 29 K is within 60 K of that of the extensively-studied 
planet HD 209458 b. 
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1 INTRODUCTION 

Of the known extrasolar planets, those which transit their par- 
ent stars comprise a rich source of information on the formation 
and evolution of gas giants. The physical properties of transiting 
extrasolar planets (TEPs) can be determined - with a little help 
from stellar e volutionary theory - to precisions appr oaching a few 
percenQ (e.g. IXorres etai]|2008l : lsouthwortlj|2009l) . The radii of 
the known TEPs are of particular interest as indicators of the in- 
ternal st ructure and heating mechanisms of irr adiated giant plan- 
ets (e.g. iFortnev et aHllOOTI : ISatygin et al]|2009i) . Also, the opaci- 
ties of different species in the atmospheres of TEPs change with 
wavelength, causing variations of their radii as measured from 
transit light curves. These wavelength-dependent variations can 
be used to probe th e atmospheric composition of the objects (see 
iBurrows et al J2007h . 



'Fortn ev et alj ^20081) divided the known gas-giant TEPs into 
two classes according to the presence of appreciable atmospheric 
opacity from oxidised titanium and vanadium. They calculated the 
observed radius of a planet with a radius of 1.20 Rjup at 1 bar pres- 
sure and a surface gravity of 15ms~^, for wavelengths of 0.3- 
30 /^m. pM-class planets have high TiO and VO opacity leading to 
temperature inversions and an increased transit-derived radius over 
a broad range of optical wavelengths. These objects are predicted 
to have apparent radii which are the smallest at the bluest optical 
wavelengths (<370 nm) and the greatest at 450-700 nm, the differ- 
ence being 4% for the 1 .20 Rjup planet and even larger for one with 
a lower surface gravity (larger atmospheric scale height). By con- 
trast, the pL planets are less highly irradiated and have negligible 
TiO and VO opacity. Their transit radii are predicted to be smaller 
and less variable, with narrow peaks around the Nal 589 nm and 
K I 768 nm doublets. A similar variation for pL planets was found 
bv lBurrows etal] jloiQ) for the case of HD 209458. 



See: |http : / /www, astro . keele . ac .uk/^ jkt/tepcat/| 
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Table 1. Log of the observations presented in tliis work. A^obs is the number of observations and 'Moon ilium.' is the fractional illumination of the Moon at 
the midpoint of the transit. The aperture sizes are the radii of the software apertures for the star, inner sky and outer sky, respectively. 



Transit 


Date 




Start time 


End time 




Exposure 


Filter 


Airmass 


Moon 


Aperture 


Scatter 








(UT) 


(UT) 




time (s) 








ilium. 


sizes (px) 


(mmag) 


1 


2010 


8 23 


20:24 


01:33 


201 


60 


Stromgren u 


1.00 -> 


2.01 


0.994 


17, 30, 50 


3.61 


1 


2010 


8 23 


20:24 


01:33 


204 


60 


Gunn g 


1.00 -s. 


2.01 


0.994 


30, 45, 70 


1.20 


1 


2010 08 23 


20:24 


01:33 


204 


60 


Gunn r 


1.00 -S> 


2.01 


0.994 


30, 45, 70 


0.93 


1 


2010 08 23 


20:24 


01:33 


199 


60 


Johnson I 


1.00 -S> 


2.01 


0.994 


27, 45, 65 


1.27 


2 


2011 05 26 


21:18 


02:47 


135 


120-100 


Gunn r 


1.38 -> 1.0 


1.07 


0.294 


16, 26, 45 


0.75 


3 


2011 07 19 


20:39 


02:04 


176 


90-33 


Gunn r 


1.03 1.0 


1.53 


0.861 


20, 30, 50 


1.53 



IZahnle et al. has shown that sulphur is an important el- 

ement for opacity in exoplanet atmospheres, and might serve as a 
probe of the metallicity of TEPs. S2 opacity becomes large at 240- 
340 nm for temperatures above 1200 K, and HS absorbs strongly 
in the wavelength interval 300-460 nm. These effects should lead 
to large apparent radii at blue-optical and ultraviolet wavelengths. 
Another relevant phenomenon is Rayleigh scattering, which is pro- 
portional to A^* and so should cause the transit radii of TEPs to be 
larger at s horter wavelengths. B oth effects run opposite to the pre- 
dictions of lportnev et al.l <2008l) : it is not yet clear which processes 
dominate the measured radii of exoplanets. 

The two well-studied TEPs, in terms of variation of radius 
over optica l wavelengths are HP 209458 b and HD 189733 b. For 
the former, iKnutson et alj j2007l) presented HST grism photome- 
try of two transits divided into ten wavelength intervals covering 
290-1030 nm. D ifferences of ra dius with wave l ength are unclear 
from these data: iBarmanI ( |2007|) andlSing et alj (l2008h found sig- 
nifican t variati ons whereas iKnutson et alj ^20071) and ISouthworthI 
( l2008h did not. ICharbonneau et alj ( I2OO2E has found evidence for 
Nal absorption in the planetary atmosphere based on an earlier 
dataset (Brown et al. 2001). These characteristics point to a pL 
classification for HD 209458 b; its moderate eq uilibrium temper- 
ature of Te'q = 1459 ± 12 K dSouthwortdboid) puts it roughly at 
the boundary between the pL and pM clas ses. For HD 1897 33 b, 
HST observations cove ring 5 50-1050 nm jPont et al.l [20081) and 
290-570 nm (Sing et alj|201ll) returned only gradual radius vari- 
ations with wavelength, which have been interpreted as indica- 
tive of Rayleigh scattering from a high-altitude atmospheric haze. 
HD 189733b has = 1191 ± 20K (Southworth 2010) so is 
firmly in the pL planetary class. There has therefore been no clear 
detection of TEP radius variation s indicative of a pM classification. 

HAT-P-5 was discovered bv lBakos et H] bOOTh and is a com- 
paratively normal TEP system consisting of a 1.2 M© star or- 
bited by a l.lMjup planet. Its planet is slightly hotter (Tcq = 
1516 ± 29 K) t han HD 209458 b s o is on the border of the pL 
and pM classes. lKrivovetalJj201lh found a flux excess at 12 and 
22 /xm indicative of the presence of dust rings at a distance fo 
several AU from the parent star. In this work we present the first 
follow-up transit photometry of HAT-P-5 obtained since the dis- 
covery paper, covering four optical passbands. We probe for radius 
variations between these passbands and improve the measurements 
of the physical properties of the system. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed one full transit of HAT-P-5 on the night of 2010/08/23, 
using the 2.2 m telescope and BUSCA imager at Calar Alto Astro- 
nomical Observatory. BUSCA uses dichroics to split the incoming 



Table 2. Excerpts of the light curve of HAT-P-5. The full dataset will be 
made available at the CDS. 



Telescope 


Filter 


BJD(TDB) 


Diff. mag. 


Uncertainty 


CAHA22 


u 


2455432.35006 


-0 


00004 


0.00316 


CAHA22 


u 


2455432.56365 


-0 


00208 


0.00582 


CAHA22 


a 


2455432.35006 





00156 


0.00119 


CAHA22 


a 


2455432.56365 





00004 


0.00135 


CAHA22 


T 


2455432.35006 


-0 


00013 


0.00094 


CAHA22 


r 


2455432.56465 





00144 


0.00099 


CAHA22 


I 


2455432.35006 





00015 


0.00129 


CAHA22 


I 


2455432.56465 





00093 


0.00134 


Cassini 


T 


2455708.38967 





00066 


0.00081 


Cassini 


r 


2455708.60776 


-0 


00120 


0.00119 


Cassini 


r 


2455761.36079 


-0 


00104 


0.00115 


Cassini 


r 


2455761.53645 





01876 


0.00195 



light into four wavelength intervals, which traverse different arms 
of the instrument and are incident onto four CCDs. This allows 
photometry to be obtained in four passbands simultaneously. The 
arms of the instrument cannot be operated individually: their focus 
and CCD settings are controlled together. The full field of view is 
12x12 arcminutes, which the CCDs sample with a plate scale of 
0.176" per pixel. 

For our observations we chose to have the dedicated BUSCA 
109-mm diameter Stromgren u and Johnson / filters, in the bluest 
and reddest arms respectively. For the intermediate arms we opted 
for 50-mm diameter Gunn g and r filters from the standard Calar 
Alto filter catalogue. The latter two filters have a better through- 
put but lead to a reduced field of view: a circle approximately 6 
arcmin in diameter. The g and r filters are also thinner than the u 
and I filters, leading to a difference in focus between the two filter 
sets. We defocussed BUSCA in such a way as to have the g and 
r observations more defocussed than those in it and /, in order to 
compensate for the lower throughput of the latter two filters. 

The amount of defocussing was set by the requirement that the 
brightest pixel within the point spread function (PSF) of the target 
and comparison stars (on all four CCDs) have not more than 35 000 
counts in one observation. An initial exposure time of 90 s was cho- 
sen, which is slightly s maller than our usual value of 120 s (see 
ISouthworth et al. l2009al ibl) to limit the impact of the bright moon. 



The CCDs were binned 2x2 in order to shorten the readout time. 
The autoguider was not defocussed along with BUSCA, so could 
be used to maintain the pointing of the telescope. 

We observed two transits of HAT-P-5 in May and July 201 1, 
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Figure 1. Plot of the residuals of the timings of mid-transit of HAT-P-5 versus a linear ephemeris. The timings obtained by amateur astronomers are plotted 
using open circles. 



using BFOSC mounted on the 1.52 m G. D. Cassini Telescop^ at 
Loiano Observatory, Italy. For more informati on on the use of this 
instru ment for defocussed photometry see Sou thworth et al. Il l20ld 
l201lh . Both transits were observed through a Gunn r filter. The 
second transit was curtailed by cloud before egress. It also suf- 
fered from the gradual drift of the telescope back into focus over 
the night, requiring the exposure time to be successively reduced 
from the original 90 s down to 33 s. As with BUSCA, the telescope 
pointing was maintained using the autoguider. A summary of the 
observational data is given in Table[T] 

One further transit was observed by NR using a 25 cm 
aperture Meade LX200 telescope sited near Sorrento, Italy (lati- 
tude 40° 37' 07.37" North, longitude 14° 21' 27.46" East, altitude 
275 m above sea level). The CCD used was an SBIG ST7 operating 
at a plate scale of 1.16" px~^, mounted behind an r filter. Data re- 
duction was performed using Maxim DL Pro ^following standard 
procedures. 

Several images were taken with the telescopes properly fo- 
cussed, to verify that there were no faint stars within the defocussed 
PSF of HAT-P-5. We find no evidence of a second star near the PSF 
of our star of interest, and conclude th at the transit is not dilu ted by 
contaminating light. This agrees with lDaemgen et al] l l2009h . who 
obtained high-resolution observations of HAT-P-5 as part of a high- 
speed 'lucky' imaging survey of twelve TEP host stars. 

Data reduction was undertaken using standard methods. Aper- 
ture photometry was per formed using the IDiI^ASTROLIeQ imple- 
mentation of DAOPHOT ( IStetsoniri987h . The apertures were placed 
by hand and shifted to account for pointing variations, which were 
measured by cross-correlating each image against a reference im- 
age. We tried a wide range of aperture sizes and retained those 
which gave photometry with the lowest scatter compared to a fit- 
ted model (see below). We find that the shape of the light curve is 
very insensitive to the choice of aperture sizes, and to whether flat 
fields are used in the reduction process. The times of observation 
were converted from UTC to barycentri c Julian date on th e TDB 
timescale, using the IDL procedures of Eastman et aljboiol) . 

Differential photom etry was obtained using a n optimal ensem- 
ble of comparison stars dSouthworth et alj|2009al) combined with a 
polynomial fit to the data outside transit. The choice of comparison 



Table 3. Times of minimum light of HAT-P-5 and their residuals versus the 
ephemeris derived in this work. 



Time of minimum 


Cycle 


Residual 


Reference 


BJD(TDB) - 2400000 


no. 


(JD) 




54241.77700 ±0.00022 


-427.0 


0.00013 


Bakos et al. (2007) 


54969.56817 ±0.001 10 


-166.0 


-0.00031 


Gregorio (AXA) 


54969.56837 ±0.00180 


-166.0 


-0.00011 


Mendez (AXA) 


54969.56947 ± 0.00090 


-166.0 


0.00099 


Naves (AXA) 


55005.81857 ±0.00080 


-153.0 


-0.00007 


Norby (AXA) 


55036.49157 ±0.00150 


-142.0 


-0.00028 


Srdoc (AXA) 


55050.43371 ±0.00118 


-137.0 


-0.00051 


Brat (TRESCA) 


55050.43453 ±0.00155 


-137.0 


0.00031 


Tmka (TRESCA) 


55334.85873 ±0.00171 


-35.0 


0.00021 


Shadick (TRESCA) 


55379.47582 ±0.00060 


-19.0 


0.00172 


Vila'agi (TRESCA) 


55387.83985 ±0.00122 


-16.0 


0.00033 


Garlitz (TRESCA) 


55432.45534 ± 0.00066 


0.0 


0.00024 


This work (u) 


55432.45436 ±0.00033 


0.0 


-0.00074 


This work (g) 


55432.45460 ± 0.00027 


0.0 


-0.00050 


This work (r) 


55432.45437 ± 0.00025 


0.0 


-0.00073 


This work (/) 


55708.5 1460 ±0.00022 


99.0 


0.00061 


This work (r) 


55482.64913 ±0.00081 


18.0 


0.00151 


Shadic (TRESCA) 


55716.87995 ±0.00086 


102.0 


0.00054 


Garlitz (TRESCA) 


55736.39824 ±0.00091 


109.0 


-0.00048 


Ayiomamitis (TRESCA) 


55744.76350 ± 0.00092 


112.0 


-0.00064 


Garlitz (TRESCA) 


55761.50126 ±0.00290 


118.0 


0.00628 


This work (NR) 



stars has a negligible impact on the shape of the observed tran- 
sits. A second-order polynomial was required for the CAHA light 
curves, but a first-order one was enough for the Cassini data. The 
resulting photometry is given in Table|2l and reaches a scatter as 
low as 0.75 mmag for the first Cassini t ransit. We al s o obta ined the 
2:-band and i?-band data presented in iBakos et al ] l l2007h for in- 
clusion in our analysis. These datasets cover two and four transits, 
respectively. 



^ Information on the 1.52m Cassini Telescope and BFOSC can be found 
at |http : / /www .bo .astro . it /loiano/] 
^ |htt p : / /www ■ cyanogen . com/| 

^ The acronym IDL stands for Interactive Data Language and is a 
trademark of ITT Visual Information Solutions. For further details see 
|http: / /www, ittvis ■ com/P roduct Services/ IDL . aspx| 
^ http://idlastro.gsfc.nasa.gov/ 



3 ANALYSIS 

We have measured the physical properties of t he HAT-P-5 s ystem 
using the methods described by ISouthworthI feOOSl l2009l (20ld 
uOIlh . thus retaining consistency with the Homogeneous Studies 
project results. A detailed description of our approach is given by 
those works, so we only summarise it briefly below. 
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Figure 2. Phased light curves of HAT-P-5 compared to the best JKTEBOP 
fits using the quadratic Hmb darkening law. The sources of the light curves 
are labelled on the plot. The residuals of the fits are plotted at the base of 
the figure, offset from zero. 



3.1 Period determination 

As a first step in our analysis we fitted all of the light curves individ- 
ually using JKTEBOP (see below). The absolute values of the obser- 
vational errors from our pipeline (which stem ultimately from the 
APER subroutine) are optimistic, so were rescaled for each dataset 
to give a reduced 'x^ ofXv = 1- The times of mid-transit were then 
measured, and their uncertainties estimated from Monte Carlo sim- 
ulations. The second Loiano transit was disc ounted from this as the 
data cover only part of the transit event (see lGibson et al.l l20091. 
To our times of minimum we added 14 timings obtained 




-0.04 -0.02 0.00 0.02 0.04 

Orbital phase 



Figure 3. As for Fig.l2] but for the z and R light curves of HAT-P-5 pre- 
sented bv lBakos et allBoOTi) . 

by amateur astronomers and made available on the AXjAQ and 
TRESC^ll 

websites. We included only timings based on data of 
low scatter and covering a full transit. All timings whose timescale 
is not stated were assumed to be on the UTC system and converted 
to TDB. 

The resulting measurements of transit midpoints were fitted 
with a straight line to obtain an orbital ephemeris. Cycle zero is 
specified to be the midpoint of the transit observed at CAHA. We 
found: 

To = BJD(TDB) 2 455 432.45510(10) + 2.78847360(52) x E 

where E is the number of orbital cycles after the reference epoch 
and quantities in brackets denote the uncertainty in the final digit 
of the preceding number. The quality of the fit is ~ 1.46, which 
implies that a linear ephemeris is not a perfect fit to the transit tim- 
ings. There is no systematic deviation from the predicted transit 
times, so we do not regard this as sufficient evidence for transit 
timing variations. The x 'i instead interpreted as an indication 
that the errorbars on the timings are too small, which we have com- 
pensated for by increasing the uncertainties in the ephemeris above. 

3.2 Light curve modelling 

The light curves were modelled using the JKTEBOljfl code as de- 
scribed in 'South worthI jlOOSt) . The main information to be gleaned 
was the fractional radii of the star and planet, ta = and 
rh = where a is the orbital semimajor axis and i?A and i?b 
are the true radii of the two objects. The fitted parameters were the 
sum and ratio of the fractional radii, rA+^b and k = and the 
orbital inclination, i. A mass ratio of 0.001 was adopted to govern 
the shapes of the biaxial ellipsoids representing the two component 
objects. Limb darkening (LD) was imposed using each of five laws 
and with three choices of whether the LD coefficients were fixed or 

^ Amateur Exoplanet Archive, http : / /brucegary . net/AXA/x . htm| 
^ The TRansiting ExoplanetS and CAndidates (TRESCA) website can be 
found at, http://var2. astro. cz/EN/ tresca/ index . php 
* JKTEBOP is written in FORTRAN?? and the source code is available at 
ihttp : / / www .astro ■ keele ■ ac ■ uk/~ jkt / 
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Table 4. Parameters of the JKTEBOP fits to the light curves of HAT-P-5. The final parameters are the weighted means of the results for the six datasets. Results 
from the literature are included at the base of the table for comparison. 



Source 


''A + ''b 




* ( ) 


''A 


'i"h 


CAHA 5-band 
v_,/a n /a i -udiiti 
CAHA 7-band 
Loiano r-band 
Bakos z-band 
Bakos i?-band 


0.1407 ±0.0080 
n 1 AAl A- ft 0077 

0.1404 ±0.0059 
0.1562 ±0.0084 
0.1464 ±0.0077 
0.1413 ±0.0095 


0.1 130 ±0.0025 
n 1 1 -k n 0099 
0.1 134 ±0.0017 
0.1 159 ±0.0026 
0.1094 ±0.0020 
0.1111 ±0.0025 


87.8 ± 1.7 

88.2 ±1.5 

86.1 ±1.1 

87.2 ±1.4 
87.4 ± 1.9 


0.1264 ±0.0069 
1 909 -t- oo^^^^ 

0.1262 ±0.0051 
0.1400 ±0.0073 
0.1320 ±0.0067 
0.1272 ±0.0084 


0.0143 ±0.0011 
01 so -1- 0010 

0.0143 ±0.0008 
0.0162 ±0.0011 
0.0144 ±0.0010 
0.0141 ±0.0012 


Final results 






87.18 ±0.61 


0.1296 ± 0.0027 


0.01467 ±0.00041 


Bakos et al. (2007) 
Torres et al. (2008) 


0.1477 
0.1478 


0.1106 ±0.0006 
0.1 106 ±0.0006 


86.75 ± 0.44 
86.75 ± 0.44 


0.133 ±0.003 
0.1333 ±0.0033 


0.01471 
0.01472 



fitted. U ncertainties were calcula ted using both Monte Carlo sim- 
ulations jSouthworth et al. | l2004 ) and a residual-permutation algo- 
rithm (as implemented bv lSouth worth 2008). The larger of the two 
values was retained fo r each output quan tity. The orbital eccentric- 
ity was fixed to zero teakos et al.ll2007h and co ntaminating light 
was taken to be negligible ( Daemgen et alj2009l) . 

The CAHA r and I data were solved individually, and the 
ti-band light curve was ignored at this point due to its much greater 
scatter. The Loiano transits were solved simultaneously but with 
the orbital period fixed at the value determined in Section lBTI We 
also considered the «-band and 7?-band follow-up data presented in 
lBakosetalJj2007h . In all six cases we found that the best solutions 
were obtained when fitting for the linear LD coefficient and fixing 
the nonlinear coefficient at a theoretically-predicted value (but per- 
turbing it by ±0. 1 on a flat distribution during the Monte Carlo and 
residual-permutation simulations). 

An overall set of photometric parameters was calculated for 
each dataset and is given in Table|4] Tables of individual fits for 
each light curve are given in the Supplementary Information which 
accompanies this work. The corresponding best fits are shown in 
Fig. [2]for the data prese nted in this work and in Fig.[3]for the data 
frorn bakos et al.l j2007h . The final photometric parameters are the 
weighted mean of the values for each dataset. At first glance the 
Loiano results are outliers, but this impression is not backed up by 
statistics: the value of the agreement of the individual parame- 
ters with respect to the weighted mean is 1.5 for k but smaller than 
0.6 for the other photometric parameters. T able[4]also sho ws a com- 

S iariso n with the results from lBakos et al.i (2007.) and Torres et al.l 
20081) . These authors agree well with our final parameter values, 
but find errorbars which are notably smaller than ours for those 
datasets in common. 



3.3 Variation of planetary radius with wavelength 

Once the final photometric parameters had been established, each 
light curve was fitted with i and ta fixed in order to investigate the 
possibility of a variation of with wavelength. The resulting rb 
values should have reliable relative uncertainties from which the 
common sources of error have been removed. We inflated the er- 
rorbars to account for systematic noise following the approach 
and w ith groups of up to ten consecutive datapoints (e.g. lWinn et al.l 
1200'^ . One issue to bear in mind when fixing i and ta is that any 
variations in these parameters, or effects arising from systematic 
noise, will be concentrated in the remaining free parameters, in- 
cluding rb. 



Whilst rb might vary detectably with wavelength, the LD 
characteristics of the parent star certainly do. In order to make any 
errors in the theoretical representation of LD by stellar model at- 
mospheres clear, we calculated rb values with the linear LD coef- 
ficient included as a fitted parameter as well as solutions in which 
both coefficient were fixed to theoretical values. An excellent dis- 
cussion of the treatment of LD in transit light curves is given by 
Howarth (2011). 

Fig-El shows the variation of rb with wavelength with the lin- 
ear LD coefficient fitted (top panel) and fixed (middle panel). The 
corresponding results are given in Table[5] Clear variations with 
wavelength are seen, and are robust against the treatment of LD. In 
all cases fixing LD coefficients leads to rb values which are smaller, 
an effect which is most pronounced for the it-band data. Changing 
the LD coefficients to their physical limits (i.e. the total LD at the 
limb of the star is zero or unity) cannot shift the it-band rb into 
alignment with the other results. 

The referee raised concerns that correlated noise in the pho- 
tometry could be influencing the observed variations in rb, and that 
by fixing i and rA to common values any systematic errors would 
be concentrated in rb. We have therefore shown in Fig.|4] (bottom 
panel) the values of rb obtained when i and ta are not constrained 
to be the same for each light curve. The differences between the 
three panels in Fig.|4] therefore highlight the effect of adopting a 
common set of system parameters. 

So is the variation of rb with wavelength real? The two r- 
band datasets are in good agreement, which is encouraging. The 
It-band result is very discrepant with those from the other bands, 
but the reality of this phenomenon is arguable. Of all optical wave- 
lengths, photometry blueward of the Balmer jump suffers the most 
from the strength and variability of Earth-atmospheric extinction, 
although the altitude of CAHA (2168 m above sea level) will mit- 
igate that to some extent, u- and f7-band photometry is also often 
handicapped by a lack of comparison stars, which is true in the cur- 
rent case (only one good comparison star compared to at least two 
for all other bands). Our CAHA observations were also taken at 
full moon (illuminated fraction 0.994) so the sky background was 
comparatively high in u and g. Fig.|2] shows that our it-band data 
have some correlated noise visible in the residuals. 

A second suggestion of correlated noise affecting the rb val- 
ues is that those for the z- and 7?-band data are significantly below 
the values found from our new data. Whilst this could be interpreted 
as indicative of temporal variation in the properties of the system, it 
is most likely to herald the presence of systematic errors in the pho- 
tometry. Further observations are needed to assess the significance 
of the rb variations we see. For the remainder of this work we pro- 
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Table 5. Fractional planetary radius values (rj, = -^Ji) found from the 
available data with the linear limb darkening (LD) coefficient either fitted or 
fixed, when the other photometric parameters are fixed at their known/final 
values. 



Source 



rb (LD fitted) 



rt (LD fixed) 



CAHA Ji-band 
CAHA g-band 
CAHA r-band 
CAHA /-band 
Loiano r-band 
Bakos z-hanA 
Bakos R-band 



0.01659 ±0.00031 
0.01483 ±0.00012 
0.01500 ±0.00009 
0.01533 ±0.00014 
0.01495 ±0.00010 
0.01434 ±0.00006 
0.01438 ±0.00008 



0.01581 ±0.00027 
0.01464 ±0.00011 
0.01484 ±0.00007 
0.01515 ±0.00012 
0.01475 ±0.00009 
0.01433 ±0.00005 
0.01430 ±0.00007 



ceed under the assumption that the variations of with wavelength 
are unconfirmed. 



3.3.1 Possible explanations of the radius variation 

If the observed variations of radius with wavelength were true, what 
would they imply? The re value from the ij-band data is substan- 
tially larger than those in the other passbands, which is in the oppo- 
site direction to the predictions of Fortney et al. (2008). Rayleigh 
scattering causes a larger radius at bluer wavelengths, which would 
fit the u result but not the more modest trend seen in t he g, r and I 
results . Similarly, the sulphur chemistry advocated bv lZahnle et al.l 
( I2OO9I) would lead to a larger radius in both u and g due to opacity 
ofHS. 

Next, one can calculate the size of the u-band radius differ- 
ence in units of the atmospheric pressure scale height, H. The rb 
at u is 7% (for fixed LD coefficients) or 1 1% (linear LD coefficient 
fitted) larger than that for the r band. This equates to a difference 
of 61 00 km or 9700 km. For HA T-P-5b 500 km (using data 

from lde Pater & Lissauerll200ll) . The radius variation i s therefore 
roughl y 12H or 19H, which is large but not excessive. I Sing et al] 
( l201ll. their fig. 14) measured radii for HD 189733 b which ap- 
proac hed 6H larger at 330 nm than that measured by Pont et all 
( I2OOS) in the 1000-1025 nm passband. Rayleigh scattering is there- 
fore a plausible expanation. 

Two further possibilities are apparent. Firstly, an unusally 
large transit depth was found for WASP- 1 2 bv lFossati etal] \20m 
and interpreted as evidence for a extended exosphere surrounding 
the planet, with the possible presence of a disc formed from ex- 
creted material. The observations have subsequently been analysed 
in the context of a bow shock arising from the motion of the planet 
through the stellar wind ( Vidotto et al. 2010; Llama et al. 2011). 
The observable consequences of this situation are that UV transits 
are deeper and begin earlier than those at redder wavelengths. A 
bow shock would be much weaker due to the larger orbital sepa- 
ration in the HAT-P-5 system (0.04079 ± 0.000 80 AU; see below) 
compared to WA SP-12 (0.02293 ± 0.00078 AUi lHebb et alj2009l : 
IChan et alj20ril) . Whilst our observations satisfy the criterion of a 
deeper transit at bluer wavelengths, they do not indicate that the u 
transit occurs earlier than the g, r and / ones (see TableO. 

The second possibility is that the stellar surface exhibited 
starspots during our observations, which are not directly detectable 
becasue they were not occulted by the planet. This would act to in- 
crease the depth of the transit. If the temperature of the spot(s) is 
close enough to the Teff of the star that they contribute significant 
flux at redder wavelengths, this would induce a wavelength depen- 
dence such that transits in the blue are deeper than those in the red. 
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0.01 65 7 
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0.01 55 r 
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0.0145 r 
0.0140 7 
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Figure 4. Variation of the fractional planetary radius (rb = with 
wavelength. Results fro m the CAHA data are shown in blue, from Loiano 
in green, and from the iBakos et alj j2007h data in red. The veilical lines 
represent the relative errors in the measurements and the horizontal lines 
show the full widths at half maximum transmission of the passbands used. 
The top panel gives the results for limb darkening fitted, the middle panel 
for LD fixed at theoretically predicted values, and the lower panel shows 
the results when the system parameters are not forced to be the same for 
different light curves. The filter designations are labelled at the top of the 
lower panel. 



Such a phenomenon would explain both a deeper M- band transit 
and a greater transit depth overall compared to those in lBakos et al.l 
However, it does not agree well with the observed trend in 
the g, r and I bands. The level of starspot activity would also have 
to be similar at both of our r-band epochs, which are separated by 
nine months, in order to explain their mutual agreement. 

The spot activity of HAT-P-5 A can be investigated using 
time-series photometry e xtending ove r a signifi cant time interval . 
The SuperWASP archive jPollacco et al. 2006; Butters et al.1l20IOl) 
lists 48 392 observations of HAT-P-5 over the years 2004 to 2010. 
We searched these data for periodicities indi cative of spot- i nduce d 
rotational modulation, using the method of iMaxted et al. I( l201lh . 
The periodograms of the data show only peaks near 1 d and 28 d, 
which we attribute to the diurnal and lunar cycles. We set a 95%- 
confidence upper limit of 0.5 mmag on rotational modulations, 
which confirms that HAT-P-5 is a quiet star. 
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Table 6. Derived physical properties of the HAT-P-5 system. For each set of physical properties the following quantities are the same: gb = 16. 7 ± 1.9 ms ^, 
PA = 0.791 ± 0.050 p0 and Tj^ = 1517 ± 29 K. 





This work 


This work 


This work 


This work 


This work 


This work 




(dEB constraint) 


t^oiaici moueis^ 


I moQeis } 


y iciaiiiu moueis^ 


V noo moQeis } 


t^LJOCr moQeis J 


K\y ( km s ) 


1 o 1 /I r» 

loi.z ih 4.y 


159.8 ±1.9 


1 1 1 o 1 

159.1 ± l.i 


158.6 ± 2.9 


icon 1 1/1 

158. U ih 2.4 


1 CO C 1 1 f\ 

158.5 ± 3.0 


Ma (Mo) 


1.240±0.113 


1.185 ±0.042 


1.170±0.051 


1.159 ±0.063 


1.147 ±0.052 


1.156±0.065 


Ra(Rq) 


1.161 ±0.044 


1.144 ±0.029 


1.1 39 ±0.029 


1.135 ±0.031 


1.132 ±0.029 


1.135 ±0.032 


logSA (cgs) 


4.402 ±0.021 


4.395 ±0.019 


4.393 ±0.019 


4.392 ± 0.020 


4.390 ± 0.020 


4.392 ±0.020 


Mb ( Mjup) 


1.10±0.13 


1.07±0.11 


1.06±0.11 


1.06±0.11 


1.05 ±0.11 


1.05±0.11 


Rh ( R jup) 


1.279 ±0.053 


1.260 ±0.039 


1.254 ±0.040 


1.250 ±0.042 


1.246 ±0.040 


1.250 ±0.042 


Pb (PJup) 


0.494 ± 0.067 


0.501 ±0.066 


0.503 ± 0.067 


0.505 ± 0.067 


0.507 ±0.067 


0.505 ±0.067 


e 


0.0580 ±0.0064 


0.0589 ±0.0063 


0.0591 ±0.0063 


0.0593 ± 0.0063 


0.0595 ± 0.0063 


0.0594 ± 0.0064 


a(AU) 


0.04166 ±0.00126 


0.04104 ±0.00049 


0.04086 ± 0.00059 


0.04073 ± 0.00074 


0.04059 ±0.00061 


0.04070 ± 0.00076 



Table 7. Final physical properties of the HAT-P-5 system, compared with results from the literature. Where two errorbars are given, the first refers to the 
statistical uncertainties and the second to the systematic en'ors. 





This work (final) 


Bakos et al. (2007) 


Torres et al. f2008) 


MaCMq) 


1.163 ± 0.065 ±0.022 


1.160 ±0.062 


+0.043 
^•^^'-0.081 


Ra(R©) 


1.137 ± 0.032 ±0.007 


1.167 ±0.049 


i.iOO_Q (,52 


logSA (cgs) 


4.392 ± 0.020 ±0.003 


4.368 ±0.028 




PA (Po) 


0.791 ± 0.050 




729+0 058 


Mb(Mjup) 


1.06 ± 0.11 ±0.01 


1.06±0.11 


1 06+"-" 

^•'-'"-0.11 


Rh (Rjup) 


1.252 ± 0.042 ±0.008 


1.257 ±0.053 


1 254+0.051 


9b (ms-2) 


16.7 ± 1.9 


16.5 ±1.9 


16.6 J g 


Pb (Pjup) 


0.504 ± 0.067 ±0.003 


0.50 ±0.08 


"sn+0 09 


Te'q (K) 


1517 ±29 




1539+32 





0.0592 ± 0.0064 ± 0.0004 




0591+° °°^'* 


a (AU) 


0.04079 ± 0.00076 ± 0.00026 


0.04075 ± 0.00076 


04071+" °°°''^ 


Age (Gyr) 


1 7+3.5 +0.4 
' -1.4 -0.6 


2.6 ± 1.8 


2 6+2-1 



In conclusion, we observe a different transit depth in u com- 
pared to g, r and I. We also fi nd lower transit dep ths for the R 
and z observations obtained by iBakos et al.l ( |2007|) . We have at- 
tempted to interpret these effects in the context of a variation of 
opacity with wavelength in the planetary atmosphere, the presence 
of a bow shock, or spot activity on the host star None of these pos- 
sibilities provide a completely satisfactory description of the data, 
but Rayleigh scattering remains plausible. Our preferred explana- 
tion is that the differences in transit depth are due to systematic 
effects in the available photometry. 



3.4 Physical properties of the HAT-P-5 system 

We determined the physical properties of the H AT-P-5 system us- 
ing the approach described in ' Southworth|j2009l) . Starting with the 
measured photometric parameters and kno wn orbital velocity am- 
phtude of the star (J^a = 138 ± 14ms~^ lB"akos et alj|2007h . we 
guessed a velo city amplitude for the planet (/S'b) and used standard 
formulae (e.g. lHilditclJllOOll) to obtain an estimate of the proper- 
ties of the system. These comprise the mass, radius, surface gravity 
and mean density for the star (Ma, Ra, log (?a and pa) and planet 
(Mb, Rh, <7b and pb), and the orbital semimajor axis (o). We then 
interpolated within the predictions of theoretical stellar models to 
find the effective temperature (Teg) and Ra corresponding to the 
estimated Ma aii d measured metallicity ([§] = 0.24 ± 0.15; 
iBakos et alj|2007h . The value of Kh was then iteratively refined 



to maximise the agreement between the estimated and observed 
Tcft (5960 ± 100 K; Bakos et al. 2007), and the estimated ^ and 
known ta. This process was performed for a grid of ages starting 
from the zero-age main sequence and stepping in 0.01 Gyr chunks 
until the star had evolved to a log 3 less than 3.5. The final result 
was a set of physical properties and a system age which gave the 
best agreement with the known Tcff and va- 

The random e rrors in this process wer e propagated via a per- 
turbation analysis jSouthworth et al.1 12005'). Our use of theoreti- 
cal models inevitably incurs a systematic error due to the depen- 
dence on stellar theory. We assessed these systematic errors from 
the interagreement of sets of physical properti es calculated usin g 
each of five independent sets of stellar models ( ISouthworthll201 ll ). 
Our final result for each physical property corresponds to the un- 
weighted mean of the five different values, whilst its statistical 
error is the largest of the five possibilities and its systematic er- 
ror is the standard deviation of the values from the five models. 
Note that three quantities are indepen dent of stellar theory s o are 
not troubled by systematic error s: Qb dSouthwort h et al.' 20071), pA 
dSeager & Mallen-Omelasll2003l) and iSouthworth 200 9). 

It is possible to avoid using stellar models by recourse to 
calibrations based on empirical measu rements of stars in eclips- 
ing binary systems jS outhworthl l2009l) . We adopted the approach 
outlined bv lEnoch et al. (2010) but w i th th e improved calibration 
coefficients calculated by Southworth (2011). The set of physical 
properties resulting from this method are shown in Table[6l along 
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with the five sets arriving from the use of stellar mod els. Our final 
proper ties are given i n Table|7] The results obtained bv lBakos et alj 
ll2007h and lTorres et al. ( 2008) are in very close agreement with our 
own, despite being based on much sparser photometric data. 



4 SUMMARY 

We present photometry of four transit events in the HAT-P-5 extra- 
solar planetary system, obtained using telescope-defocussing tech- 
niques and reaching scatters as low as 0.75 mmag per point. One of 
these transits was observed in four passbands simultaneously, using 
the BUSCA imager on the CAHA 2.2 m telescope. We used these 
data to improve the measured orbital ephemerides and physical 
properties of the system. HAT-P-5 is well-characterised, although 
it would benefit from further observations to refine the spectro- 
scopic orbit and atmospheric properties of the host star. The planet 
is slight ly too large to match the predicted radii of coreless gaseous 
planets jFortnev et alj2007l : lBaraffe et al.l2003) . confirming its sta- 
tus as an inflated hot Jupiter. HAT-P-5 becomes the sixtieth TEP 
system to be stud ied as part of the Homgeneous Studies project 
jSou thworth 2008) and its physical properties have been placed in 
the Transiting Extrasolar Planets Catalogu^fl 

Both components of HAT-P-5 have similar temperatures to 
those in the prototype transiting system HD 209458, so the planets 
are expected to have similar atmospheric characteristics. We used 
our multi-band photometry to search for variations of the measured 
planetary radius with wavelength, which could indicate the pres- 
ence of Rayleigh scattering or the opacity of specific molecules. 
We find that the radius in the u band is significantly larger than 
at the other optical wavelengths, by either 12 or 19 pressure scale 
heights depending on the treatment of limb darkening when fitting 
the light curves. We also find that the planet radius measured in 
earlier R and 2-band data is smaller than that from our data. These 
phenomena are most likely due to the presence of systematic errors 
in the photometric data. Alternative possibilities include Rayleigh 
scattering and temporal changes in the system properties, but ex- 
planations involving starspots or a bow shock do not match existing 
observations. 
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